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The polymerization dynamics of double-stranded polymers, such as actin fila- 
ments, is investigated theoretically using simple chemical kinetic models that ex- 
plicitly take into account some microscopic details of the polymer structure and the 
lateral interactions between the protofilaments. By considering all possible molecular 
configurations, the exact analytical expressions for the growth velocity and disper- 
sion for two-stranded polymers are obtained in the case of the growing at only one 
end, and for the growth from both polymer ends. Exact theoretical calculations 
are compared with the predictions of approximate multi-layer models that consider 
only a finite number of the most relevant polymer configurations. Our theoretical 
approach is applied to analyze the experimental data on the growth and fluctuations 
dynamics of individual single actin filaments. 
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I. INTRODUCTION 

Cytoskeletal proteins such as actin filaments, intermediate filaments and microtubules 
are rigid multifilament polymers that play a variety of roles in biological systems, includ- 
ing organization of cell structures, transport of organelles and vesicles, cell motility and 
repro duct ion. Biological functions of these proteins are mostly determined by the pro- 
cesses that take place during their polymerization. However, our understanding of the 
coupling between these biopolymer's structure and functions and their growth processes is 
still very limited. 

In recent years the number of experimental investigations of the growth mechanisms and 
dynamic properties of rigid multifilament biopolymers at a single-molecule level increased 
significantly.— '^'S'i'^'Siifliii Dynamic behavior of individual microtubules have been character- 
ized by a variety of experimental techniques such as video and electron microscopy, fluores- 
cence spectroscopy, and optical trap spectrometry ,^i^i^i'^i^i^ whereas the studies of the single 
actin filaments have just begun . Many unusual phenomena in the assembly dynamics 
of these biopolymers have been observed, such as treadmilling for microtubules and actin 
filaments,'^'i2iii and microtubule dynamic instability.— 

Recent experimental investigations of the single actin filament length fiuctuationsi^ raised 
many questions about the actin polymerization dynamics. A large discrepancy in kinetic 
rate constants estimated by average length change in the initial polymerization phase and 
from the analysis of length fluctuations in the steady-state phase (a factor of 40) has been 
observed. One of the possible explanations for this discrepancy might be the oversimpli- 
fied theoretical model, used in analysis,-^ that neglected the polymer structure and lateral 
monomer-monomer interactions at the actin filament tips. Similar problems have been found 
in the growth dynamics of individual microtubules under the influence of external forces . ^1^1^^ 
For example, the depolymerization rate constants determined from a phenomenological de- 
scription of microtubule dynamics, that do not incorporate the biopolymer's structure and 
lateral interactions between the filaments, contradict to values measured in some indepen- 
dent bulk chemical kinetic experiments.-^ 

A large volume of experimental results stimulated many theoretical investigations of 
polymerization dynamics for rigid multifilament proteins. In one approach, it was sug- 
gested that the growth of rigid biopolymers is controlled by thermal fluctuations.-^^*^^^ This 
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is a basic idea of the so-called polymerization ratchet models. In a different approach, 
more phenomenological chemical kinetic (stochastic) models have been used to describe the 
biopolymer's growth dynamics , ^^i^' i 'i^^i^^ In these simplified phenomenological models it is as- 
sumed that the overall dynamics is a balance between polymerization and depolymerization 
processes, however the microscopic details of the polymer structure and the differences in lat- 
eral inter-subunit interactions are not taken into consideration. In addition, the microtubule 
assembly dynamics has been studied extensively by computer simulations . ^^i^^i^^i^^ 

Recently, we introduced a set of simple stochastic models for the description of the growth 
of rigid biopolymers consisting of protofilaments, that explicitly includes the geometric 
structure and monomer- monomer lateral interactions.^^ According to this approach, only 
few polymer configurations are relevant for the growth dynamics because of inhomogeneity 
in lateral interactions that comes from the geometry of polymer ends. These configurations 
are selected using the criteria that the distances between the protofilament tips in each 
configuration should be less than nd, where n = 1, 2, ■ ■ ■ and is a polymer subunit length. 
For n = 1, in the so-called "one- layer" approximation, the mean growth velocity and a 
dispersion (or a diffusion coefficient) of polymer's length have been calculated exactly for 
any number of protofilaments and for the arbitrary geometry of biopolymer's growing end. 
The quality of the "one-layer" approximation has been discussed for the simple case of the 
growth of polymer made of two protofilaments, i.e., N = 2, which is closely related to 
actin filaments. In this case, a full dynamic description that includes all possible polymer 
configurations provided exact expressions for the mean growth velocity. It was found that 
the predictions of the approximate model are approaching the exact results for large but 
realistic values of the lateral interactions. 

Although the approach presented in Ref.— provided a very good description of exper- 
imental results of individual microtubule growth under external forces, and it suggested 
a reasonable way of coupling of the microscopic structure of the biopolymers with their 
dynamic properties, there are several theoretical problems with this method. First, the 
criteria that helps to determine the finite number of the most relevant configurations in 
"n-layer" approximations is introduced using an arbitrary cut-off distance, but it would be 
more desirable to derive it from more fundamental thermodynamic and kinetic arguments. 
Second, the method used for calculating exactly the full dynamic model of the growth of 
two-stranded polymers allowed only the determination of the mean growth velocity, and not 
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the dispersion. However, the simultaneous knowledge of the dispersion and the velocity is 
crucial for understanding the growth mechanisms of rigid biopolymers.— The goal of this 
paper is to address these theoretical issues by studying the growth dynamics of two-stranded 
biopolymers with more detailed microscopic description. In addition, we will analyze the 
experimental observations on the growth of single actin filaments. 

This article is organized as follows. The dynamics of the two-stranded polymers that 
can grow only from one end is discussed in Section II, while the description of the polymer 
growth from both ends is given in Section III. The application of the developed chemical 
kinetic models for the experiments on single actin filaments is presented in Section IV. The 
results are discussed and summarized in Section V. The mathematical details of calculations 
of dynamic properties are given in Appendix. 

II. GROWTH DYNAMICS OF ATTACHED TWO-STRANDED POLYMERS 

Consider a growing rigid polymer that consists of two protofilaments. The polymer is 
attached to a nucleating site on a surface and may polymerize or depolymerize only from 
the free end, as shown in Fig. 1. This is the attached two-stranded polymer.— The building 
block of this polymer is a monomer subunit of length d. The lateral offset between two 
parallel protofilaments is equal to a and it can generally vary between and d. For actin 
filaments the monomer size is d = 5.4 nm and the lateral shift is a = d/2 = 2.7 nm.^^*^ 

There are infinite number of possible polymer configurations that differ from each other 
by the geometry of the growing end and the total length of the polymer. For labeling these 
configurations we use a pair of integers that count the number of monomer subunits in each 
protofilament. Without loss a generality, let us choose a first (bottom) protofilament as a 
leading one in a configuration (/, /), where the distance between protofilament tips equal to 
a: see Fig. 1. The second (top) protofilament will be the leading one in the configuration 
labeled as {1,1 + 1), where the distance between the tips is equal to d — a. Generally, as 
illustrated in Fig. 1, for configurations (/, l — k), where A; = 0, 1, 2, ■ ■ ■ , the distances between 
the parallel protofilaments are equal to a + kd and the first protofilament is the leading one, 
while for configurations (/ — A;, / + 1) the distances between protofilament tips are d — a + kd 
and the second protofilament is the leading one. 

All dynamic transitions in the attached two-stranded polymers may be characterized by 
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a set of four pairs of transition rates that depend on the local geometry of the growing end. 
As shown in Fig. 1, the attachment and detachment rates from the leading protofilament 
for all polymer configurations are given by uq and Wq, respectively. Transition rates ui and 
Wi describe the polymerization and depolymerization events when the overall length of the 
polymer docs not change. For example, as presented in Fig. 1, the transitions between the 
configurations {1+1,1) and {1 + 1,1 + 1) are given by these rates, while in both configurations 
the polymer length is d{l+l)+a. A monomer can attach with the rate us to the protofilament 
2 of the configuration {1,1), or it can detach with the rate ws from the same protofilament 
of the configuration {1,1 + 1). Similarly, the monomer subunit can associate with the rate 
Ui^s to the protofilament 1 of the configuration (/,/ + 1), or it can dissociate with the rate 
Wis from the same protofilament of the configuration (/ + !,/ + 1): see Fig. 1. Note that 
the subscript indexes (0, 1, S and 1 — S, where S — a/d) are equal to the fractions of the 
lateral bond between the monomers on the parallel protofilaments created or broken in a 
given transition. 

The overall kinetic scheme of the system, that includes all possible states and transitions, 
is shown in Fig. 2a. The polymer growth dynamics can be described by a set of master 
equations for each configuration. Solutions of these equations, that characterize the steady- 
state growth dynamics of an attached two-stranded polymer, are outlined in a full detail in 
the Appendix A. Below we present only the final results for the mean growth velocity and 
dispersion. 

Within the full dynamic description, the growth velocity appears formally as a sum of 
two terms, namely. 



V = Vo{l-(3)+Vi 



(1) 



where 



U5U1S - WsWis 



(2) 



Us + Wi + Ui_s + Wis 



and 



Vi = d{uQ - wqP). 



(3) 



The parameter /3 (0 < /3 < 1) is given by 

^ _ Uq + Wi 

Ui + Wo 

The expression for the diffusion coefficient is also consists of two terms. 



(4) 



D^Do{l-p)+D 



(5) 
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where the terms are given by the following expressions, 

(f {usui_s + wswi_s - 2Al) 



and 



2 US + WS + ui-s + Wis 

2{Ao + wo){uo - Ao(3) 



Ui + Wq 



The auxiliary function Aq is defined as 



Ao = Vo/d. 



(6) 



(7) 



(8) 



The dynamic properties of the growing polymer strongly depend on the lateral interac- 
tions between the parallel protofilaments. It can be seen from the fact that the transition 
rates for binding or unbinding the monomer subunit are directly related to a lateral inter- 
action free energy per monomer Qh via the detail balance conditions. It can be shown thafeSl 



Us/ws = ^^o/wo7^^ s = 0,5,1- 5,1, (9) 

where 

'J = exp{-gh/kBT). (10) 

This observation is the basis for the approximate theoretical description of the growth 
of rigid biopolymers with filaments, the one-layer model, that we developed earlier.— In 
this model, only polymer configurations with the distances between the protofilament tips 
less than the monomer length d are considered. It allows then to calculate the mean growth 
velocity and dispersion explicitly. Specifically, for = 2 it can be shown that 

Vone-layer = aud Done-layer = Dq. (H) 

Since the full dynamic description of the growth of the attached two-stranded polymer is 
now available, the quality of the one-layer approximation can be easily checked in this case. 
As shown in Fig. 3, the one- layer model describes the growth dynamics reasonably well for 
large lateral interactions. In the limit of infinite lateral interactions the predictions from the 
approximate theory become exact. However, the convergence of the approximate one-layer 
results to exact quantities strongly depends on the geometry of the growing polymer end, 
specifically, on the ratio between the lateral shift a and the monomer length d. The best 
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description can be obtained for the symmetric case 5 = a/d = 1/2, while for other geometries 
this approximation is less successful. 

One of the advantages of the one-layer approximation is its ability to be easily extended 
to include more polymer configurations. As a better approximation, it is natural to con- 
sider configurations where the distances between the protofilament tips do not exceed 2d, 
i.e., a two-layer model. The number of polymer configurations is still finite and the dy- 
namic properties can be easily calculated following the approach presented for the one-layer 
approximation.— In the two-layer approximation, for two-stranded polymers it can be shown 
that 

and 

Dtwo-layer = ^ _^ ^ + -^l)) (13) 



where ^ ^ 

(14) 



2 



, „ 2{A + wo)iuo-A(3) 
Uo + WoP 



Ui + Wo 

with the parameter A given by 

A = Vtwo-layer / d. (1^) 

The results of the two-layer approximation for the mean growth velocity and dispersion are 
also presented in Fig 3. The agreement with exact full dynamic properties is very good 
even for weak lateral interactions [gn > 2kBT), and the results do not depend much on the 
specific geometry of the growing polymer end. 

The multi-layer approach can be used to describe the growth dynamics of any rigid poly- 
mer consisting of parallel protofilaments. It is important to understand thermodynamic 
and kinetic justifications for this approximations. The fact that the growth dynamics for 
two-stranded polymers (A^ = 2) can be analyzed exactly is very useful to make the connec- 
tion with n-layer approximate description. The kinetic schemes for full dynamic description 
and for the one-layer and two-layer approximations for two-stranded polymers are shown 
in Fig. 2. Comparing different kinetic diagrams, we can see that the one-layer approach 
corresponds to the main chemical pathway, while the two-layer approximation also takes 
into account the closest branched states. Thus, the ra-layer approximations can be thought 
of as a series expansion (with n terms) of full dynamic description, where the value of energy 
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of lateral interactions determines how good is the expansion. The higher the lateral inter- 
actions the smaller number of terms is needed in order to describe successfully the growth 
dynamics of rigid multi-filament polymers. 

III. GROWTH DYNAMICS OF FREE TWO-STRANDED POLYMERS 

Now consider a two-stranded polymer that can freely grow from both ends. Define Xiit) 
and XFj{t) as time-dependent coordinates of the "left" and "right" ends of the polymer. The 
growth velocity is defined as 

Vfree = J^{[xRit)-XLm, (16) 

and the angular brackets mean averaging over all possible growth pathways. It can be easily 
seen that the mean growth velocity of a free polymer is a difference between two terms, 

Vfree = J^i^t)) " ^(^^W) = Vr ' (17) 

where Vr and Vl are one-end growth velocities, determined explicitly in Eqs. dH), 0, and 
®. 

Similarly, the expression for the dispersion is given by 

Dfree = ~ [{[xR{t) - X.it)]') - {[xn{t) - X,{t)]f] . (18) 

It can be simplified into the following equation, 

Dfree = ~ [{xRitf) + {xLitf) - {xnit))^ - {xL{t)f - 2(x^(t)Xi(t)) + 2{x R^t)) {x Lit))] . 

(19) 

The polymerization dynamics at both ends are independent from each other, that means 
that 

{xR{t)xL{t)) = {xR{t)){xL{t)). (20) 

This leads to the conclusion that the dispersion of free growing polymer can be presented 
as a sum of two one-end dispersion terms, 

Dfree = DR + DL, (21) 

where the explicit expressions for Dr and are given by Eqs. (0), (0) and ((7j). 
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IV. APPLICATION OF CHEMICAL KINETIC MODELS FOR THE 
DESCRIPTION OF EXPERIMENTS ON INDIVIDUAL ACTIN FILAMENTS 

The growth dynamics of the single actin filaments has been studied experimentally using 
the fluorescence microscopy with total internal reflection.i^ The assembly of actin filaments 
was observed during the "polymerization" phase, i.e. at initial stages of the process, as well 
as the steady-state conditions, when the growing rate of the barbed end was compensated 
by the shortening at the pointed end. The main finding was that the kinetic rate constants 
estimated from the length change in the single actin filaments for the initial period differ 
considerably from that estimated using the length fluctuation analysis in the steady-state 
phase. The set of rate constants measured at the initial poljTiierization conditions mainly 
agrees with values obtained in other experimental studies,— i^Si^ while the rate constants 
estimated from the measurements of length fluctuations at steady-state conditions were 
30-45 times higher. 

Several possible explanations for this discrepancy has been suggested.-^ The flrst one is 
that the kinetic constants obtained in the steady-state phase may be intrinsically different 
from those obtained in the initial phase of polymerization. The change of nucleotide com- 
position of the growing, or shortening tips in the time course of the process was indicated 
as a probable cause. The depolymerization velocity of ADP-actin is known to be an order 
of magnitude higher than that of ATP-actin at the barbed end. The second possible reason 
for the discrepancy in the rate constants may be due to the possibility that the "effective" 
size of polymerization-depolymerization unit may not necessarily correspond to a monomer. 
The authors speculate that one plausible way to eliminate this divergence is to set "effec- 
tive" size of unit 5-6 times higher. However, this contradicts to widely accepted picture that 
the elementary step in the growth of actin fllaments is adding or removing a single actin 
monomer 

In analyzing the experimental data on growth dynamics of single actin fllaments the 
simplifled phenomenological picture has been used.-^ Here, we investigate another possibility 
to explain the difference in the kinetic rate constants by using a chemical kinetic model with 
better description of polymer ends geometry and chemical interactions between monomers. 

To estimate the parameters that describe the growth dynamics of actin fllaments we note 
that (5 = a/(i = 1/2, and, using the detailed balance conditions [see Eqs. the rate 
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constants can be written in the following form 

Us = ui_s = ^^o7^'/'"*"^''^ Ui = Uo'yf^+'^, 

ws = Wis = ^07'^^''^""^^^, wi = wo'y-^^'^. (22) 

Coefficients fi/2 and /i reflect the different values of activation energies for specific polymer- 
ization and depolymerization events. Although the exact values of these parameters cannot 
be measured experimentally, they might be estimated quite reahstically as —0.5 < /1/2 < 0.5 
and — 1 < /i < 1. It implies that the sub unit attaches faster to the site where the stronger 
lateral contact is created. Similarly, the detachment is slower if a stronger lateral bond 
should be broken. For simplicity, in our calculations we consider only the case 

/1/2 = /i = 0, (23) 

and, as we checked, for other values of these parameters the results do not deviate much 
from the one presented here. 

Equations ()22|) imply that the growth dynamics of actin filaments can be described by 
using only 3 parameters: Uq, Wq and 7. The parameters Uq = k^C (where C is the concen- 
tration of free actin monomers in the solution) and Wq are the association and dissociation 
rates from the leading protofilaments, and 7 is a measure of lateral interactions in actin 
filaments. The mean growth velocity for each end of actin filaments can be presented in a 
simple form, 

V=^{uo-Woh){i'/'-i-'/' + 2), (24) 
while the dispersion of the polymer length at each end is given by more complex expression. 



D 



2 



]{uo + woh){i'l' + 4 - 57-^/^ + 2/7) - (1 - 2r"' + 1/7)^^^?^ 
4 M0 + W0/7 



(25) 



In order to apply our explicit expressions to describe the single actin filaments growth 
the elongation rate constants for each end of the polymer should be known. However in 
the single-molecule experiments by Fujiwara et al. W the growth dynamics of each end sepa- 
rately has not been measured. Nevertheless, for calculations we can use the data from other 
investigations where the polymerization dynamics at both ends has been characterized quan- 
titatively in the similar experimental conditions.— Kinetics of actin polymerization for the 
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barbed end can be described as V'^ = k'^C — with fc^ = 11.6 /iM~^s~^ and = 1.4 
s~^. Comparing this phenomenological expression with the exact one [see Eq. (j211)] allows 
us to estimate the parameters uq and wq. Here we also use the realistic estimate of energy 
of lateral interactions (7/1 ~ 6 ksT^"^ that gives 7 ~ 400. These parameters are then applied 
to compute the contribution to the dispersion from the barbed end using Eq. Similar 
approach is utilized for the pointed end, for which the mean growth velocity can be described 
phenomenologically as V'^ = k^C — k^_ with fc^J. = 1.3 /iM~^s~^ and k^_ = 0.8 s~^.^^ As a 
result, the overall dispersion of single actin filaments at steady-state concentration Cq = 0.17 
fiM can be estimated as -D ~ 1.0 x 10^'^ yum^/min. It should be noted that this procedure 
depends weakly on the value of 7. Also in these calculations we used the subunit length 
d = 5.4 nm, and lateral off-set a = 2.7 nm. 

In the single actin filaments experimentsi^ the measurement of fluctuations at steady state 
conditions produced the dispersion of D ^ 1.1 x 10~^ yum^/min, which is approximately 10 
times larger than the value calculated above. The difference is significant and it implies 
that the chemical kinetic models with detailed description of polymer ends and monomer- 
monomer interactions still cannot explain fully the experimentally observed fluctuations in 
growing actin filaments. However, our theoretical treatment does not take into account the 
hydrolysis of actin- ATP monomers and related processes. It might be expected that these 
processes can significantly effect the growth dynamics of actin filaments. 

It is interesting to note that the dispersion for actin filament assembly (with a = d/2) 
is a non-linear function of monomer's concentration, as shown in Fig. 4. At high concen- 
trations of actin monomers the dispersion is proportional to concentration, while for low 
concentrations there is a weak deviation from linearity. This dependence contrasts to the 
observed and calculated behavior of the mean growth velocity. It will be interesting to 
measure experimentally the concentration dependence of dispersion. 

V. SUMMARY AND CONCLUSIONS 

We investigated theoretically the growth dynamics of two-stranded polymers where as- 
sociation and dissociation of monomers can take place from both ends. Because the poly- 
merization events at each end are independent from each other, we argued that the overall 
polymer elongation dynamics can be described as a combination of growth processes at each 
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end separately. 

For attached rigid two-stranded polymers, that made of two protofilaments and can only 
elongate from one end, we developed a chemical kinetic model of the growth. This model 
takes into account the exact relative positions of two protofilaments and both lateral and lon- 
gitudinal chemical interactions between the monomers. Considering full dynamic chemical 
kinetic scheme, the exact and explicit expressions for the mean growth velocity and dis- 
persion have been derived in terms of rate constants of binding and unbinding of monomer 
subunits. Because of the geometry of the polymer end and the monomer-monomer inter- 
actions, the growth properties of two-stranded polymers depend only on three parameters, 
namely, the rate constants of attaching or dissociating from the leading protofilament and 
the energy of lateral interactions. 

The exact full dynamic description of the growth of two-stranded polymers, that accounts 
for all possible configurations, has been compared with a set of n-layer approximate models 
that consider only the most relevant polymer configurations. It was shown that the approx- 
imate approach is successful because it captures the main features of full dynamic kinetic 
diagram. In addition, the approximate description becomes better for larger lateral inter- 
actions between the monomer subunits. It has been concluded that n-layer approximations 
might be viewed as a series expansion of the full dynamic description of polymer growth 
dynamics. It implies that the approximate approach can be used to describe the growth 
dynamics of rigid biopolymers with many protofilaments, like microtubules or intermediate 
filaments. 

The full dynamic chemical kinetic model of the growth of two-stranded polymers has been 
applied to analyze the experimental observations on single actin filaments growth. Using the 
kinetic rate parameters and the realistic estimate of the lateral interactions extracted from 
bulk chemical kinetic measurements of actin filaments, we calculated the overall dispersion 
in the length fluctuations of single actin filaments. The obtained value of the dispersion 
was approximately 10 times smaller than the experimentally observed.-^ The difference is 
significant and it implies that other processes, not accounted by current theoretical analysis, 
contribute to the dispersion of the single actin filaments. It was argued that this discrepancy 
is due to the fact that the hydrolysis in the polymer molecule is not accounted in our 
theoretical approach. 

In addition, we also discussed the concentration dependence of dispersion. Our theoret- 
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ical calculations suggest that the dispersion of actin filaments depends weakly non-linearly 
at low concentrations of free monomers, and it approaches the linear dependence at large 
concentrations. It will be very important to measure the concentration dependence experi- 
mentally since it will give a valuable information on the mechanisms of growth and it will 
provide a direct check of the validity of our theoretical picture. 

In a future, we plan to investigate the effect of hydrolysis of the monomers, associated 
with ATP or its analogs, on the growth of biopolymers consisting of N parallel rigid protofil- 
aments. For actin filaments {N — 2) it seems reasonable to extend the current chemical 
kinetic model, however for biopolymers with larger number of protofilaments, such as mi- 
crotubules and intermediate filaments, the coupling of hydrolysis with n-layer approximate 
approach, probably, is the most realistic approach. 
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Appendix A: Calculations for full dynamic model of the growth of attached 

two-stranded polymers 

Let us introduce the probabilities P(Z, I — k;t) and P{1 — k,l + l;t) of finding the two- 
stranded polymer in the configurations {1,1 — k) and {I — k,l + l) respectively at time t. Here 
l,k = 0, 1, ... and the two parameters in brackets correspond to the number of subunits in 
the first and second protofilaments respectively. These probabilities (at k = 0) satisfy the 
following master equations. 
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dP{l,l;t) 
Jt 



ui-sP{l - 1, 1; t) + wsP{l, l+l;t)+ uiP{l, l-l;t)+ wqP{1 + 1, /; t) 
-{us + wi_s + Uo + Wi)P{l, I; t), 



(Al) 



dP{l,l+l;t) 
Jt 



usP{l, I; t) + wi^sP{l + l,l + l;t) + uiP{l - 1, Z + 1; + wqP{1, l + 2;t) 
-{ui^s + W5 + ui + 'Wo)P{l,l+'i-;t). (A2) 
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These equations describe a set of special states on the main pathway in a chemical kinetic 
scheme, see Fig. 2a. For A; 7^ we have 

dP{l,l-k-t) ^ _ij_k.t\^ y^^pn^ / + 1 _ fc; t) + uiPil, l-l-k;t) 

at 

+WoP{l +l,l-k;t)- (uq + W0 + U1+ Wi)P{l, l-k;t), (A3) 



and 

dP{l-k,l + l]t) 



UoP{l -k,l)+ WiP{l + 1 - A;, / + 1; t) + UiP{l - 1 - fc, / + 1; t) 



dt 

+WqP{1 -k,l + 2;t)- {uq + W0 + U1 + Wi)P{l -k,l + l;t). (A4) 

The conservation of probabihty leads to 

+00 / +00 +00 \ 

J2 E l~k;t) + J2Pil-kJ+l;t)] =1, (A5) 

1=0 \k=0 k=0 / 

at all times. 

Following the idea of Derrida,— we define four sets of auxiliary functions {k = 0,1,...), 

+00 

5mW = I]^(^,^-^;^) (A6) 



1=0 

+00 



+00 

Bk,i{t) = Y,Pi^-k,l + l]t) (A8) 

1=0 

+00 

Ck,i{t) = 5Z + l)P{l -k,l + l;t) (A9) 
1=0 

where 6 = a/d. Note that the conservation of probability gives us 

+00 1 

J]J]5fc,(t) = l. (AlO) 

fc=0 i=0 

Then from master equations ()A1|) . ()A2|) . ()A3|) and ()A4|) we derive for k = 

dBofiit) ^^^^^ ^ ws)Bo^i{t) + {ui + wo)Bifi{t) - {us + wi^s + uq + wi)5o,o(t), 

(us + Wi^s)Bo,oit) + (ui + Wo)Bi^i{t) - + WS + UO + Wi)5o,i(t),(All) 



dt 
dBo,iit) 



dt 

while for k ^ {i = 0,1) it is given by 
dBf; (t) 

— — = (mo + wi)Bk-i,i{t) + {ui + wo)Bk+i,i{t) - {uo + W0 + U1 + wi)Bk,i{t). (A12) 
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Similar arguments can be used to describe functions Ckflit) and Cfc,i(t). Specifically, for 
A; = we obtain 

^'^"'"^^^ = {ui-5 + ws)Co,i{t) + {ui + Wo)Ci,o(^) - {us + Wi_s + Uo + Wi)Co,o{t) 

+ [6u,_s - (1 - 6)ws]Bo,,it) - woB.^oit), (A13) 

= {us + wi_s)Co,o{t) + (ui + Wo)Ci^i{t) - {ui_s + WS + UO + Wi)Co,i{t) (A14) 

at 

+[(1 - 5)us - Swi.s]Bo,oit) - woBi^i{t). (A15) 

For k ^ (i = 0, 1) the expressions are 

^'^''^y^ = {uo + Wi)Ck-i,i{t) + (mi + Wo)Cfc+i,i(t) - (mo + Wo + Ml + 
or 

+Mo5fc_i,i(t) - u»o5fe+i,*(t)- (A16) 

Again following the Derrida's method,— we introduce the ansatz that should be valid at 
large times t, namely, 

BkM bk,i, CkAi) ^ cik,it + Tfc,, (z = 0, 1). (A17) 
At steady state dBk,i{t)/dt = and Eqs. (IMT|l and (1X121) yield for A; = 

= + ws)bo^i + (mi + wo)bi^o - {us + wi^s + uq + wi)bofl, 
= {us + Wi_s)bo,o + (ui + Wo)^i,i - (mi-5 + WS + UO + Wi)&o,i) (A18) 
while for /c 7^ (i = 0, 1) we obtain 

= (mo + + (ui + wo)bk+i,i - (uo + W0 + U1 + wi)bk,i. (A19) 

The solutions of Eqs. (|XT8| and (IM9|l can be written in the following form 

bk,i = for {i = 0, 1), (A20) 

1 + a 

where A; = 0, 1, .., and 

Ui-s + Ws Ui+ Wo 

To determine the coefficients a^^j and j from Eq. ()A17|) . the ansatz for the functions 
Ck,i is substituted into the asymptotic expressions ()A13jl . ()A14j) and ()A16jl . yielding for 

k = 0, 

= {ui^s + ws)ao^i + {ui + wo)aifi - {us + wis + mo + m^i)«o,o, 

= {us + Wi_s)aofi + {ui + Wo)ai^i - {ui_s + ws + Uq + Wi)ao,i. (A22) 
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At the same time, for A; 7^ (2 = 0, 1) we obtain 

= (mo + wi)ak-i,i + (mi + Wo)ak+i,i - (mq + Wq + Ui + ^1)0^,^. (A23) 

The coefficients j satisfy the following equations (for A; = 0) 

ao,o = (mi-<5 + ws)To^i + (mi + ^0)^1,0 - {us + wis + uq + wi)Tofi 

+ [6ui-s - (1 - S)w5]bo,i - woh^o, (A24) 

and 

ao,i = {us + Wi-s)Tofl + {ui + Wo)Ti^i - {ui_s + ws + Uq + Wi)To,i 

+ [(1 - 5)us - 5wi^s]bo,o - wo&i,i- (A25) 

For /c 7^ (i = 0, 1) we have 

ak,i = {uo + wi)Tk-i,i + {ui + wo)Tk+iA - {uq + wo + ui + wi)Tk,i + uo&/fc-i,i - wo&fc+i,i- (A26) 

Comparing Eqs. ()A18|l and ()A19|1 with expressions ()A22jl and ()A23j) . we conclude that 

ak,i = Ahk,r, (« = 0,1), (A27) 

with the constant A. This constant can be calculated by summing over the left and right 
sides in Eq. ()A27|) and recalling the normalization condition ()A10|) . The summation over 
all ttk^i in Eqs. flA22|l and ()A23jl produces 

+00 1 

A = = [{l~5)us-Swi^s]bo,o+[Sui^s-{^-S)ws]bo^i+Uo-Wo{l-bo^o-bo^i). (A28) 

k=0 i=0 

Thus we have A = Aq{1 — (3) + Ai , where 

Ao = [{1 - 6)us - Sw^^s]bo.o + [Su.^s - (1 - 5)H&o,i = "^'"l^'' , (A29) 

and 

A^=Uo- Wo(3. (A30) 

Note that Aq does not depend explicitly on 5. 

To determine the coefficients T^, ^ , we define for all k the following function 

n = Tk,o + Tfc,i; (A31) 
flfe = ak,o + akX, (A32) 
bk = bk,o + bk,i. (A33) 
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Then we define 

Vk = (ui + wo)Tk - {uo + Wi)Tk-i. (A34) 
Now Eqs. ()A24j) . ()A25jl and ()A26j) can be rewritten as 

yo - y-i = ao + wobi - Ao{l - (3), (A35) 
Vk - Vk-i = flfc - uoh-i + wobk+i, (A36) 

with y_i = and = 1, 2, .... The solutions for these equations are given by 



Vk 



Ao{l - (3) h, - 1 j + u^h (A37) 



Summing up o or i separately for all fc, we obtain the relationship between the param- 
eters To^i and Tq q, i.e., 

6o,o[(l-5)^^5-5w^i-5] -«^o(l-/5)/(l + a) 
-fo,! — tt-fo.o H ; ) (^Ad«j 

and for sum of all we have 

J:T, = To + ^^^^^. (A39) 

It is now possible to calculate explicitly the mean growth velocity, V, and dispersion, D, 
at steady-state conditions. The average length of the polymer is given by 



\k=0 1=0 k=0 1=0 / 

+ 00 

= dJ2[CkAt)+Ck,im- (A40) 



k=0 

Then, using Eq. ()A27|) . we obtain for the velocity 

'+00 +00 

^ l(i) >= dA i 

t^oo dt 



V = \\m ^ < l{t) >=dA\y2 bk,o + ^-^.1 ) = (^41) 

\fc=0 k=0 J 

A similar approach can be used to derive the expression for dispersion. We start from 

(+00 +00 +00 +00 \ 

+ '^)'^(^' + + ^)'^(^ ~k,l + l;t)] . (A42) 

A;=0 Z=0 k=0 1=0 / 
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Then, using master equations ()A1|) . ()A2|) . ()A3|) and ()A4|) . it can be shown that 

hm 4: < l^it) > = - ^)u5 - 6wi-s]Co,o + [Suis - (1 - 5)ws]Co,i 

+ [(1 - 6yus + 6^wi^s]bo,o + [S^ui^s + (1 - Sfws]bo,i + uq 

+00 +00 

+m;o(1 - 6o,o - &o,i) + «o E [Ck,o+Ck,i] + E [Cfc,o+C^fc,i] f • (A43) 

A:=0 k=l J 

Also, the following expression can be derived using Eq. ()A40|) 

, +00 

Inn - (< m >2) = 2rfM^ [Cfc,o + Ct,i]. (A44) 

k=0 

The formal expression for dispersion is given by 

D = l lim I (< /(t)^ > - < l{t) >') . (A45) 

Then substituting into this expression Eqs. ()A43|) and ()A44|) we obtain 

D = ^ {[(1 - - 5w^-s]To,o + [5u^-5 - (1 - 5)H^o,i + - <5)2m5 + 5^wi^s]bo,o 
+l[5'^ui_5 + (1 - 5)2t(75]6o,i + luo + ^wo{l - 6o,o - bo,i) + wo(^o,o + ^o,i) 

+ 00 ^ 

+ {uo -Wo-A)Y: [Tk,o+Tk,i] > . (A46) 

fe=0 J 

Finally, after some algebraic transformations, we derive the final expression for the disper- 
sion, D = Dq + Di, which is given in Eqs. ©, (jH)) and ((Tj) in Sec. II. Note that a constant 
To cancels out in the final equation. 
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Figure Captions: 

Fig. 1 Different configurations of tfie growing two-stranded polymer molecule. The 
molecule is attached at the left end and it can grow only from the right end. The size 
of the monomer subunit is d, while a is a shift between the parallel protofilaments. The 
rates and labels for different configurations are explained in the text. 

Fig. 2 Chemical kinetic schemes for models of the growth of two-stranded attached poly- 
mers: a) full dynamic description; b) two-layer approximate model; c) one-layer approximate 
model. 

Fig. 3 Comparison of the exact dynamic properties of polymer growth calculated in the 
full dynamic description with the approximate results from one-layer and two-layer models 

for different geometries, a) Ratio of exact and approximate mean growth velocities as a 
function of lateral interactions, b) Ratio of exact and approximate dispersions as a function 
of lateral interactions. 

Fig. 4 Dispersion as a function of free actin monomers for growing single actin filaments. 
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Figure 3b 
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